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Abstract. Using the quantum Hamiltonian for a gravitational system with 
boundary, we find the partition function and derive the resulting thermody- 
namics. The Hamiltonian is the boundary term required by functional differ- 
entiability of the action for Lorentzian general relativity. In this model, states 
of quantum geometry are represented by spin networks. We show that the 
statistical mechanics of the model reduces to that of a simple non-interacting 
gas of particles with spin. Using both canonical and grand canonical descrip- 
tions, we investigate two temperature regimes determined by the fundamental 
constant in the theory, m. In the high temperature limit (kT 2> m), the 
model is thermodynamically stable. For low temperatures (kT <C m) and for 
macroscopic areas of the bounding surface, the entropy is proportional to area 
(with logarithmic correction), providing a simple derivation of the Bekenstein- 
Hawking result. By comparing our results to known semiclassical relations we 
are able to fix the fundamental scale m. Also in the low temperature, macro- 
scopic limit, the quantum geometry on the boundary forms a 'condensate' in 
the lowest energy level (j = 1/2). 



1. Introduction 

The statistical mechanics of the gravitational field will provide answers to some 
of the most intriguing questions on the relation between quantum and gravitational 
physics. These questions center on key issues such as the dominance of the semi- 
classical states, the states of geometry on the smallest possible scales (including 
microstates of black holes), the final stages of stellar collapse, and the early dy- 
namics of the universe. While general gravitational statistical mechanics remains 
relatively unexplored, considerable progress has been made in understanding quan- 
tum black holes. 

In the context of string theory, many properties of extremal and near extremal 
black holes have been derived. (See Ref. [Q for a review.) Within the context of 
four dimensional quantum gravity there has also been some progress. These have 
included studies of canonical ensembles for vacuum spacetimes in which one fixes 
the average mass of the system and computes thermodynamic quantities using the 
partition function. For Schwarzschild black holes, such canonical ensembles were 
considered in, for instance, Kastrup Q and Makela and Repo (who also consider 
the statistical mechanics of the exterior region) ||. The partition function in both 
these cases is divergent. In statistical mechanics, one can introduce a chemical 
potential to tame the divergence of the canonical ensemble. This was done by Gour 
for black holes with area spectra given by integer multiples of a fundamental area 
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The present work makes use of canonical quantum gravity using real connection 
variables ||, or, succinctly, "spin net gravity"^] Spin networks, discovered in this 
context by Rovelli and Smolin || , are linear combinations of Wilson loops based 
on Ashtekar's new variables Spin networks are one of the central elements of 
this approach to the quantization of gravity. There is a primer on spin networks 
in Ref. || and a review of the general approach in Ref. |J. The mathematics has 
evolved considerably and now provides rigorous underpinnings for the quantization 
(lO) - |lj| of the kinematic states of quantum gravity. 

Perhaps the most striking result of this held is the new picture it provides of 
geometry on the fundamental scale. In this picture, 3-geometry is represented by 
a graph. A state of quantum geometry is built from holonomies of the gravita- 
tional connection along a collection of edges. These edges are "field lines of area" 
quite similar to electric field lines in electrostatics. This new "quantum geome- 
try" is vastly different from our garden-variety geometry of meter sticks, particle 
accelerators, and regions neighboring black holes. Length frig! , area ( |]^, poj ) , 
volume ( ) , and angle |26) have been found to have discrete spectra. In 
this approach to quantum gravity, geometry is fundamentally discrete. 

Within the framework of spin net gravity, there have been a number of entropy 
calculations (|27j -Q) which are largely based on the "area ensemble". The idea is 
that one uses the area, instead of energy, as the fixed macroscopic parameter. The 
resulting ensemble is defined with respect to the area and its conjugate intensive 
quantity, the "area temperature" or surface pressure. 

We take a different approach. Gravitational systems defined in bounded spatial 
regions generically have Hamiltonians associated with their boundaries. We for- 
mulate gravitational statistical mechanics using the quantum Hamiltonian and the 
usual partition function: a sum over all states of the Boltzmann weighting factor 
e~@ H . This is only possible when the spatial manifold has boundary. In this case, 
the Hamiltonian is the quasilocal energy associated with the boundary. While our 
statistical mechanical treatment is familiar, the statistical sum is in stark contrast 
to familiar textbook presentations of statistical mechanics, where quantized matter 
fields are compared to classical geometry via the infinite volume limit. The ensem- 
bles we describe contain no matter nor is there a background geometry to compare 
to; this is the statistical mechanics of fundamental geometry. 

Even without background geometry, we nevertheless can investigate three, sep- 
arate limits to the theory: high and low temperatures and the macroscopic limit. 
The temperature regimes are determined by the scale of the theory, which is roughly 
the Planck scale. To investigate spatial boundaries with macroscopic areas, it is 
necessary to study the model when there are a very large number of intersections 
between the boundary and the underlying spin network graph. Combining both 
this macroscopic limit and the low temperature regime, we find that the entropy is 
proportional to area (with logarithmic correction) and that the system condenses 
in the lowest energy level. 

These results come from a confluence of a number of features of this model. First, 
we make use of a quantum Hamiltonian of the system |3]J| . It has a very simple 
action on spin network states. Second, spin net gravity effectively reduces the 
study of quantum gravitational field theory to quantum gravitational mechanics. 



x By spin net gravity we mean that the discrete structure of spin networks not only describes 
the states of 3-geometry but also supports dynamics in Lorcntzian four dimensional spacetime. 
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Our model is equivalent to a system of geometric particles. Third, there are clearly 
defined scales in which the theory must match previous results. This enables us to 
fix the fundamental scale of the theory. 

In the next section, we define the model that forms the basis for our study. We 
give some of the necessary classical theory before reviewing the quantum theory and 
spin network states. At each stage we give the context to the assumptions behind 



the model. These assumptions are also collected in the final subsection (2.3). In 
Section |3|, we study the canonical ensemble in the hot flat space limit and define 
the grand canonical ensemble in Section || This allows the study of the model in 
low temperature and macroscopic limits. Finally, the work is summarized and some 
implications are discussed in Section |5|. 

2. Defining the Model 

We begin the study of general gravitational statistical mechanics with a simple 
model. To motivate the assumptions and to provide the necessary background, we 
discuss the classical theory, including boundary conditions, the quantum environ- 
ment, the quantization of the Hamiltonian, and the scale of the theory. A summary 



of the model is given in 2.3 



2.1. Classical Hamiltonian Theory. We consider 4-dimensional spacetimes of 
the form M = SxK, where the spatial region S has boundary. While there may be 
more than one component of the boundary in general, such as inner and asymptotic 
components, we focus entirely on a single, 2-sphere interior boundary, <9£. We do 
not specify whether this is an intersection of the space with a horizon. 

In the Hamiltonian formalism the classical fields on the phase space of general 
relativity can be expressed using connection variables || [?J . These connections are 
the configuration variables. The real, su(2)-valued connection is A l a (x) where a 
is a spatial index and i is the index associated to the internal su(2) space. The 
conjugate momenta, or "electric field" E al (x), is a geometric variable related to 
the inverse spatial metric q ab by E al E b% — \q a b\q ah . The determinant of the spatial 
metric is denoted \q a b\- The kinematic phase space of quantum gravity is similar to 
SU(2) Yang-Mills theory. The key difference is the absence of background geometry. 

Set in a spatially bounded region, classical Hamiltonian gravity provides a clear 
selection criteria for the quasilocal energy: The (3 + l)-action must be functionally 
diffcrcntiablc. This ensures that histories remain within the phase space. In general, 
the method of functional differentiability generates boundary conditions on the 
phase space variables, gauge parameters, and surface terms [|2], |33| .[] Depending 
on the nature of the boundary conditions, surface terms may be added to the action 
to make the theory well-defined. These surface terms, without which the theory 
would be inconsistent, are the fundamental observables. 

The Hamiltonian of the theory is the surface observable arising from Hamiltonian 
constraint and is enforced by the lapse function J\f. This surface observable is the 
quasilocal energy and reduces to the ADM energy in asymptopia J3^] . It is also the 
surface term which satisfies the correct algebra: The surface term generated by the 
Hamiltonian constraint satisfies the same algebra as the constraint itself 



2 There is some confusion in the literature on this point. While functional differentiability can 
be satisfied when the lapse Af tends to zero on the boundary, this is not a necessary condition. 
Vanishing lapse on the boundary is not the only way in which to ensure functional differentiability 
of the action. 
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For these reasons, we use the surface term arising from the Hamiltonian constraint 
as the definition of the quasilocal energy and call it the Hamiltonian. 

In the real connection variables, the variation of the connection in the Hamilton- 
ian constraint generates a surface term which must either vanish when boundary 
conditions are applied or be canceled by the surface term 

Hot, (W) = [ d 2 xe l ^ k Un a AlE b ^E ak . (1) 

The lapse M has a density weight of — 1. When the lapse is non- vanishing on BY,, 
the surface term of Eq. (|l|) must be added to the Hamiltonian constraint. This 
classical expression reduces to the quasilocal energy of Brown and York Q] , to the 
ADM energy in asymptopia, and to the Misner- Sharp mass in spherical symmetry 
(32) . This Hamiltonian is the surface term which ensures functional differentiability 
of Lorentzian general relativity in terms of real connection variables. 

With the addition of the Hamiltonian of Eq. (|l|) the boundary conditions re- 
quired by functional differentiability of the Hamiltonian constraint are simply ex- 
pressed. Let n abl — Me l ^ k E a ^E hk , then 5n a K ahl \d-£ = 0. However, a complete set 
of boundary conditions includes kinematic boundary conditions, which arise from 
functional differentiability of the Gauss, diffeomorphism, and Hamiltonian con- 
straints, and dynamic boundary conditions. Dynamic boundary conditions ensure 
that the boundary conditions are preserved under evolution. These are found by 
computing the constraint algebra. See Ref. |3^] (or, for a more general analysis, 



Ref. |33f) for details. 

One set of kinematic boundary conditions which ensures functional differentia- 
bility of the constraints with the surface term of Eq. (Q) is 

S(n a e ijk E a ^E bk )\ as = 

<W|as = 0, the lapse is fixed, and (2) 
N a \gYi — 0, so that spatial diffeomorphisms vanish on the boundary. 

To ensure that the theory is dynamically well-defined one must impose further 
boundary conditions on <9E. One consistent set is 

D b K aH \ e v = 

A*|aE = (3) 
<9aA/]<3£ = so that the lapse is fixed to be a constant on dT,. 

We choose the constant lapse to be 1. 

Both sets of boundary conditions, Eqs. (||) and (||), comprise a complete, con- 
sistent set of boundary conditions for the model. As is clear from these conditions, 
the model is defined for hypersurface orthogonal spatial boundaries with fixed elec- 
tric field. Finally, we note that these boundary conditions are not in the category 
of boundary conditions which includes isolated horizons. Isolated horizons require 
that the lapse vanish on the horizon |3(i| |. 

2.2. Quantum Domain. To investigate quantum statistical mechanics, one must 
investigate both the Hilbert space of the bounded theory and the action of the 
quantum Hamiltonian. We give an outline of the state space framework here though 
the reader is encouraged to refer to references ]To|-[T^| for the elegant theory behind 
this sketch. The quantum degrees of freedom are captured by holonomies along 
paths embedded in S. Gravitational analogs of Wilson loops, a holonomy of the 
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connection gives an element of SU(2). This simple observation is the key to the 
quantum configuration space. One defines a generalized connection to be a map 
from each path to the group. Like the holonomy of the classical connection, the 
generalized connection describes parallel transport along paths. It is convenient to 
describe a collection of paths by a set of edges of a single graph. With the definition 
of generalized connections and a graph, one can introduce a measure on the space 
of generalized connections simply by using the ST/ (2) Haar measure on every edge. 
It turns out that the space of generalized connections is the completion of the space 
of smooth connections E3,0. This space is the quantum configuration space. A 
basis on this space is described by spin networks p2| . 

An open spin network s consists of the triple (G, n, i). The graph G is embedded 
in the 3-dimensional space and is comprised of a set of edges and vertices. Each 
edge is assigned a nontrivial, irreducible representation of SU(2). On every edge e 
this representation is labeled by an integer n e from the set n where n — 2j. The 
minimum representation, or edge label, is n = 1. The edges intersect at vertices. 
The number of edges incident to a vertex is called the valence. In this model 
valence-one vertices, the "open edges," only occur at the intersection of the graph 
with the boundary. Intertwiners, i, assign to each vertex v a vector i v in the tensor 
product of the representations labeling the edges incident to the vertex v. If a 
vertex is in the interior, the vector i v is invariant under the action of SU(2). If 
the vertex lies on the boundary, i v is not invariant under SU(2). We label these 
boundary intertwiner vectors by an integer which corresponds to the m,j values 
familiar from angular momentum. 

In more picturesque language, edge labels record the number of spin-1/2 "threads'' 
in a single edge of a spin network. Intertwiners label the ways in which all incident 
threads are connected. Vertices at the ends of the open edges are labeled by the 
state of the spin-1/2 threads. 

Notice that any open spin network determines a set of valence-one vertices, which 
determine a Hilbert space given by a product of angular momentum states | j rrij). 
These are the states on the boundary of the space that determine the quasilocal 
energy. 

We base the statistical mechanics on the quantization of the boundary Hamil- 
tonian of Eq. (Q). Using loop techniques and recoupling theory, the boundary 
Hamiltonian of Eq. (Q) can be quantized j3lj . The quantization of the Hamilton- 
ian is not free of choices. However, the most general form of the operator is, for 
transversal edges, 

Hqt. I a) = m J2 \/n v {n v + 2) \ s). (4) 
■uedsnG 

This operator depends only on vertices in the boundary. The lapse Af v is the scalar 
(vanishing density weight) lapse at the vertex v. In light of the boundary conditions 
of Eq. (H) , the lapse is constant on the boundary. We set J\f v = 1 for all boundary 
vertices v. 

There are two properties of the Hamiltonian operator of Eq. (Q) which impact 
the statistical mechanics. First, each edge contributes to the energy independently 
of the other edges. This simplifies the model considerably, as the total quasilocal 
energy is the sum of the contributions of each edge. Thus, as we will see, the 
problem reduces to a gas of nearly identical particles. Second, as is clear from the 
form of the Hamiltonian operator of Eq. (H), the energy of the system does not 
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depend on the vectors of the boundary vertices; the rrij value does not count. Thus, 
the system is degenerate. 

There is an important caveat to this degeneracy. We define the model only 
through the component of the system which determines the energy, the boundary 
theory. If one were to include considerations of the interior state, one would include 
the vast degeneracy, presumably infinite, of the interior intcrtwiner bases, knot 
states of the edges, and other diffcomorphism invariant aspects of the spin network 
graph @. 

In statistical mechanics, the macroscopic behavior depends critically on the par- 
ticle statistics. In our model the intersections of the surface with the network 
define a set of independent, distinguishable particles with spin. We take particles 
with different spin to be distinguishable. To motivate this it is best to describe 
measurements on the fundamental geometry. 

In quantum geometry, each intersecting edge determines the 'local' shape of the 
surface. Not only does a label n determine the area contribution of a single edge, 
but the vertex also can determine the local curvature, as in Ref. J3(J. As observed 
from the interior, the labels on these edges - the amount of geometric flux - are 
measurable via the area operator. Fluctuations in these labels are, in principle, 
observed as localized sources of radiation. The expectation is that, even if the 
macroscopic geometry is symmetric, the quantum geometry will be fluctuating on 
the microscopic scale. The localization can be achieved with the spin network itself 
if the spin network is sufficiently complex. Thus, simply based on the abstract 
graph, its labels and its intertwiners, we expect that it is possible to distinguish 
spin network edges and their fluctuations. This result may also be derived by 
examining the role of diffeomorphisms on the boundary ffij] . 

With these considerations the following picture emerges: Gravitational statistical 
mechanics of the model is based on a set of non-interacting particles. Two geometric 
particles are distinguishable if they have different spins. The particles induced from 
an edge with label n have a degeneracy of n + 1. 

2.3. Scale, Parameters, and Summary. The constant m in Eq. (|^) sets the 
fundamental scale of the theory. It is the same scale as that of the length £ which 
appears in the geometric observables. For instance, the area contribution of a single 
edge e is given by £ 2 y/n e (n e + 2). The two parameters m and I are inversely related 
with m = H/£. 

We retain the overall scale parameter in terms of m or £. There are two reasons 
for doing this. First, we wish to derive the relation between the fundamental 
scale of the model and the Planck length lp := y/ HG / c 3 . Second, this enables us 
to address a subtlety arising in the definition of the connection. As Immirzi has 
emphasized, in the canonical transformation used to define the connection variables 
there is a family of choices generated by one non-zero, real parameter 7, 1 A % a — 
r a - 7-fQ, y E ai = (l/-f)E al We absorb 7 into the constants m and L Finally, 
we retain Boltzmann's constant k throughout. 

Since m sets the scale, the regimes of interest are determined relative to this 
constant. The significant dimensionless parameter is the ratio of this energy scale 
to the thermal energy, (3m. Temperatures of order unity (T ~ 1/ km) are roughly 
10 32 K, so garden-variety energy scales of the universe are far below the natural 
scale of the theory; for quantum geometry, the universe is a chilly place. Thus 
the more realistic regime is the low temperature limit. Throughout this paper we 
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also use such terms as "microscopic" and "macroscopic." By microscopic we mean 
boundaries with area on the order of £ 2 . By macroscopic we mean boundaries with 
areas much greater than I 2 , i.e. boundaries on scales of elementary particles up 
through astrophysical objects. As we see in Section ||, to use microscopic results in 
the macroscopic regime requires a study of the scaling of physical quantities as they 
diverge. All semiclassical results need to be closely examined in low temperature, 
macroscopic limits. 

When we compare the thermodynamics of this model with semiclassical results, 



we fix the fundamental scale of the theory. This is done in Section 4.2, where 
we compare the area-entropy relationship of our model to the Bckcnstcin-Hawking 
area-entropy relation. 

In summary, the remainder of the paper deals exclusively with the model defined 
by: 

1. A spatial manifold having boundary <9E. 

2. We assume that the quantum boundary conditions are a straightforward im- 
plementation of the classical boundary conditions of Eqs. || and ||. 

3. A quantum Hamiltonian given by Eq. ^ with: 

(a) Fixed, unit lapse on dT,, 

(b) A degeneracy of n e + 1 associated to every edge e. 

4. We assume valence-one boundary vertices, i.e. the intersections G n <9£ occur 
along single edges 

5. We assume edges with different spins can be distinguished. 

Hence, the model reduces to a non-interacting, degenerate gas of distinguishable 
particles with spin. 



3. Canonical ensemble 

This section provides groundwork for both the canonical and grand canonical 
ensembles. In the canonical ensemble, we concentrate on the high temperature 
limits of the theory. 

Given any boundary with non-vanishing area, a set of edges will intersect the 
boundary. In the canonical ensemble, we fix the total number of intersections 
between the spin network graph and the boundary to be N, while allowing the 
energy of the system (the labels on the edges) to fluctuate. These intersecting edges 
are the only edges in the spin network that are open. 

We investigate the statistical mechanics of N distinguishable non-interacting 
geometric particles. To specify a state of the system, we need two things: a list 
of edge labels {n^} for those edges that intersect the boundary <9£, and an "to/' 
value for each of these edges. Since the particles are distinguishable, {n.i} will be 
an ordered list. 

The canonical partition function, Z = Tr e~@ H , becomes 

Z N {T) = Y,9{{ni})e-^ + -+^) ( 5 ) 

with the quantum Hamiltonian of Eq. (Jij) with unit lapse. The energy contribution 
of the zth edge is given by e ni := m$/rii{ni + 2), where to is the fundamental mass 
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scale for the system. The degeneracy due to the boundary vertices is 

N 

9({ni}) = Y[(ni + 1). (6) 
As the particles are independent, the partition function reduces to 

>+l)e-^ . (7) 

With this partition function it is an easy matter to compute the thermodynamic 
potential 

J*f> (8) 

where 

oo 

/:=t> + l)e-*» (9) 



is the one particle partition function. 

Normally at this stage one would study the behavior of the partition function 
and thermodynamic quantities in the infinite volume limit. As we are studying 
the underlying geometry itself, we cannot do this. We have no background for 
comparison. Nonetheless, we can explore the high temperature limit. 

In the high temperature limit, (5m <S! 1, the sum in / may be approximated by 
an integral^] which can be integrated exactly 

/oo 

2 4 . _ (10) 

= 777^7 e ~ VI/3m ( 6 + 6\/3(3m + sV3(/3m) 2 + V¥{l3mf 
(pro) 4 V 

The single particle partition function behaves as 

12 3 12 

{(imp 2 (pm) 

for high temperatures. In this temperature regime, physical quantities are easy to 
compute. The entropy is 

BA 

S = = kN (In 12 - 41n(/3m) + 4) . (12) 

The heat capacity simply scales with the number of intersections 

G = T df = ^ (13) 
As C > 0, Le Chatelier's principle is satisfied and the system is thermally stable. 
One initially would expect that since the behavior of Hawking radiation requires 
that the specific heat of black holes to be negative, this high temperature result 
would preclude black holes. However, there are quite general arguments [B7| that 
the specific heat of Lorentzian black holes ought to be positive. Further, York 
found that, using Hawking's result for spherically symmetric black holes, black 



f^T^-^ + OWrn)*-— (11) 



3 More precisely, one can establish upper and lower bounds on the integral using the Cauchy- 
Maclaurin integrals. In the high temperature limit, both these bounds have the leading order 



behavior displayed in Eq. 
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holes in a box can be thermodynamically stable. In fact, York found that there were 
two possible mass solutions, with one stable and one unstable. He described these 
as two different phases of the system, one with stable black holes and one of hot 
flat space. To map out the correspondence between the model York considers and 
the present model, we would need to add an additional boundary and a parameter 
which measured the size of the outer boundary. The topology of our present model 
prevents us from seeing both solutions. 

Average energy and area are also easy to compute. The average energy is given 

by 

N 

(E) 

i=l 

with 



0e r , 

f ,' 



Since all the expectation values for each of the particles are all identical, (E) = 
N(Ei) for any i. In the high temperature regime we make the approximation 

/>oo 

(Ek) « ^ / dx(x + 1) ^x(x + 2)e' 0m ^+ T \ 
J Ji 

The integral can be done exactly. The leading order behavior gives a total average 
energy of 

(E) = ANkT (14) 

in the high temperature limit. Similarly, the area of the bounding 2-surface, in the 
high temperature limit, is given by 

(A) = 20£ 2 N f^j . (15) 

All of these physical quantities are proportional to the number of intersections 
N. From the entropy ( |l2| ) and the area ([if]) results, we see that the entropy 
and area obey the Bekenstein-Hawking proportionality S — k-^ only at specific 
temperatures. In this high temperature approximation, this occurs when (3m « 
0.78 and 5.8. Since the Bekenstein-Hawking result is valid in lower temperatures 
(J3m 1), these are spurious solutions. 

It is perhaps not surprising that semi-classical physics differs with these results 
at high temperature scales. We should not be surprised to discover that the semi- 
classical calculations, derived in the low temperature domain, do not agree with 
results derived for temperatures much greater than the Planck scale! 

As is often the case in statistical mechanics, to investigate lower temperatures it 
proves easier to work with the grand canonical ensemble. This is what we turn to 
next. 

4. Grand Canonical Ensemble 

In the statistical mechanics of particles the grand canonical ensemble includes 
varying particle number. In this case we allow the number of intersections between 
the spin network state and the boundary to fluctuate. This means that we consider 
the category of dynamics in which spin network states evolve by changing the 
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underlying graph in a manner consistent with kinematical constraints and boundary 
conditions.^ 

A new element in our approach is the use of the chemical potential. Though 
used in the Schwarzschild case in Ref. B, the chemical potential has not been 
used within the context of spin net gravity. Allowing the number of particles to 
fluctuate provides a number of benefits. The scale of the theory, to, is so large that, 
in the canonical ensemble, any change of the system is fantastically unlikely for 
temperatures such as occur in our universe. The dynamics appears to be "frozen" . 
Indeed, the analysis would seem to require that all geometric objects have only 
two spin-1/2 edges! The chemical potential allows interesting behavior to occur 
at lower temperatures and for macroscopic configurations. Further, the chemical 
potential provides a measure of the geometric interaction between the interior and 
the boundary. 

Calling the number of intersections, A, and the chemical potential, the re- 
sulting grand canonical partition function is 



Z(7» 9({ni})e-^ + - + ^-» N ), (16) 

N=l {m} 

where, as before, g({rii}) is the degeneracy of the edge states. Since we consider 
boundaries with non- vanishing area, the sum over A begins at 1^| Each edge 
independently contributes to energy, so the partition function simplifies. 

2 = E (E( n+1 ) e ^ (en " M ) 




l-fz 

in which the fugacity z is defined as usual, z = e^, and / is defined in Eq. (||). 

This partition function diverges as fz approaches 1. There is physics in this 
divergence. To investigate the divergence we first compute some physical quantities. 

The thermodynamic potential is given by 

n=-kT(lafz-]n(l-fz)). (18) 
The ensemble average of the number of intersections is 

For physical configurations, the ensemble average (A) must be greater than or equal 
to 1. The lower bound, (A) = 1, introduces a relation between temperature and 
chemical potential, 

In/ 

» = - — 

which selects physical configurations when the boundary has non-vanishing area. 
A plot of the relation (A) = 1 is given in Fig. (Q) . 

4 Note that the classical boundary conditions, Eqs. (g) and ^ allow the boundary metric, 
including "area density" n a E at , to change. 

5 Readers familiar with spin networks will see that a single edge intersection would violate the 
admissibility conditions, in particular the evenness condition. One can account for this in the 
ensemble without affecting the physical results presented here Jisj |. 
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Figure 1. The restriction that the average number of intersec- 
tions be greater or equal to 1, (N) > 1, yields a relation between 
the temperature and the chemical potential. The upper curve is 
the plot of the relation, (N) = 1. There is a maximum tempera- 
ture for every chemical potential. The lower curve, or "macroscopic 
curve," is the relation fz = 1 which corresponds to infinite inter- 
section number. The curve near zero temperature is roughly linear, 
with a slope of about —0.7 and an intercept of ~ \/3- Both curves 
have an error on the order of 10~ 3 at the Planck temperature 
(3m = 1. 

The number of intersections diverges when the product fz approaches 1. Since 
/ is a monotonically increasing function of the temperature, these restrictions yield 
a maximum temperature for every chemical potential. This relation is the upper 
curve in Fig. (Q) . Temperatures higher than shown in the figure send the partition 
function to infinity and (N) to unphysical, negative values. All physical configura- 
tions lie inside the region between the two curves in Fig. ([!]). We call this region 
bounded by the curves (N) = 1 and fz ~ 1 the "macroscopic band." 

To better understand the geometry at the divergence, when fz — ► 1, it is worth 
deriving the area of the bounding surface. The area may be obtained through the 
simple relation between the Hamiltonian and the area at each intersection i plfl 




(20) 



A short calculation shows that 



( A ) = zH N(fz) N - + l)e"- 



N=l J 71=1 
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Thus, 

(A)={N){a), (21) 

in which (a) is the single particle average area. 

Since the area scales with (A), macroscopic areas occur for large (N). The 
divergence (A) — -> oo when fz — > 1 is physical in that the divergence selects the 
region of parameter space where macroscopic objects live. In this model of discrete 
space, macroscopic objects are found in (l//3m, n/rn) space between the curves 
shown in Fig. ([!]). 

At low temperatures, it is easy to characterize the macroscopic limit. When 
(3m 3> 1, fz may be expressed as 

fz = 2e-^-"> ^1 + J2 (^) e _/3(en " ei) ^ = 2e-^-^ + 0(e~^ mA ), (22) 

where A := (y8 — v 7 ^) is the numerical difference in the first two energy levels. 
Neglecting the higher order terms, the macroscopic condition fz ~ 1 gives us 

m 2 

M--— (23) 

This relation, the low temperature, leading order approximation to the top curve 
in Fig. ([!]), proves useful when comparing the entropy and area. We continue to 
use the "~" whenever we express quantities in the low temperature, macroscopic 
limit. 

4.1. Condensation. By examining the occupation numbers of the edge labels, one 
can find macroscopic occupation of the lowest level. If one rewrites the partition 
function in terms of number of intersections at the fcth level, A&, then one quickly 
finds the average 

_ 1 dlnZ _ (fc + l) e -fr fc 

As might be expected, these quantities diverge with (A) for all k. It is more 
interesting to investigate how the relative proportions (Nk)/(N) behave. The ratio 
is 

M = *+l e -** (24 ) 
(AO / K ' 

which, curiously, is independent of the fugacity z. With the same techniques just 
employed in the limit for fz in Eq. (|2^) , at low temperatures the occupation of the 
lowest possible level, k = 1 behaves as 

||Ul-0( e -<^). (25) 

As (N) becomes macroscopic, so does (N\). The system "condenses" in the spin- 
1/2 state, meaning that in the limit of large (N) and at low temperatures, the 
occupation of the lowest state accounts for virtually all the particles. This is much 
like Bose-Einstein condensation, which is characterized by a finite fraction of states 
in the zero momentum state in the infinite particle limit 

lim M > o. 
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If the ratio vanishes in the limit, then there is no condensation. While the analogous 
ratio in the quantum gravity case, Eq. (|2^), is positive and near unity, it is difficult 
to be precise as the (N) — > oo condition is the macroscopic limit and results in a 
more complicated relation between /i and (3. (The condition fz ~ 1 ties the chemical 
potential to the temperature, producing the macroscopic curve.) Nevertheless, it 
is easy to characterize the system in the low temperature limit; the model relaxes 
exponentially into the lowest state as the ratio of neighboring low lying states 

(Ni) 2 

vanishes exponentially as T — ► 0. 

There is another difference between the condensation in this model and in a 
Bose-Einstein system. The fluctuations, (AN 2 ) = (N 2 ) — (N) 2 , in intersection 
number behave as 

<AA 2 ) = fz 



{l-fz) 2 ' 

Hence, the relative dispersion (AN 2 ) / (N) 2 goes to I in the low temperature, macr- 
coscopic limit. The relative dispersion of the number of intersections at the lowest 
level also behaves this way. In a Bose-Einstein system, the analogous relative dis- 
persions would vanish at low temperatures. In our model, despite the fact that the 
spin-1/2 state is macroscopically occupied, we still have a tempest of fluctuations. 
As discussed in Section || these fluctuations lead to potentially observable effects. 

4.2. Gravitational thermodynamics. To check whether this model of quantum 
gravitational statistical mechanics reproduces the results of gravitational thermo- 
dynamics, we first derive the remaining macroscopic parameters of the system, 
entropy and average energy. Since we have not specified the boundary conditions 
so carefully as to pick out a black hole solution, it is interesting to see whether 
the relations of black hole mechanics are included in our thermodynamics. The 
remaining physical quantities are derived in a straightforward way, using standard 
statistical mechanics. 

Defining the ensemble average of energy of a single intersection as 

oo 

— ^e„(n+l)e-^, 



the energy is given by 



The entropy is given by 



(E) = _^|p) =(w) e (26) 



S= "(lf) =<:1 °W + *ta(M + ^W-fW. (27) 

The first thermodynamic relation we study is the Bekenstein-Hawking formula 
S = kA/Alp. A quick glance at the equations for area, entropy, and energy in this 
model show a more complicated dependence than this simple proportionality. How- 
ever, this is to be expected as the thermodynamic result holds at low temperatures 
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and macroscopic geometries. At low temperatures, the area of Eq. ( pl| ) is simply 
(N) times the basic area element, 

(A) ~ £ 2 V3{N), 

so that the area scales with the number of intersections. 

Using the relation between chemical potential and (3 at low temperatures, Eq. 
(p3|), we find that 



V3 t 2 V p 

in the low temperature, macroscopic limit. This establishes the proportionality 
between entropy and area, with a logarithmic correction. This is identical to the 
result for black holes! 

The reason for the proportionality factor ln2/-\/3 is easy to see. Suppose the 
surface area is quantized with uniformly spaced levels A = nod 2 , where n = 1,2,..., 
a is a fixed parameter, and I is the fundamental length. This suggests that the 
bounding surface is like a quilt sewn together of Planck scale patches, each of area 
al 2 . If every patch has two possible states, then a surface of area A = nal 2 has 
2™ surface states. This degeneracy contributes to the entropy and S = In 2" = 
hi2(A/al 2 ). This is precisely what we have in Eq. (pif), with a = s/3. Thus, the 
entropy is equal to the log of the number of states of the lowest level times the 
number of units of the fundamental area element \fil 2 . 

In the context of non-perturbative quantum gravity, a logarithmic correction 
has been derived before by Kaul and Majumdar [Q. However, they found that 
the correction term was negative, indicating that the semiclassical result was the 
maximal entropy. In other derivations of the logarithmic term, the corrections arise 
from quantum effects involving matter fields on a classical background. 

To reproduce the thermodynamic result of Bekenstein, we fix the fundamental 
scale of the theory. We find that 

e = lM'- (29) 

Matching with the thermodynamic result also fixes our fundamental energy scale.[] 
It is also interesting to investigate the mass-temperature relation for spherical 
black holes to see whether these solutions also are found in this class of gravitational 
statistical mechanics models. The expected relation is given by 

(M) oc i 

Since the energy is simply proportional to the intersection number and the temper- 
ature, it seems that this model appears not to include black hole solutions of this 
type. 

While the results for spherically symmetric black holes do not match the ther- 
modynamics of this ensemble, the entropy bound Jil| 

- < 2ttRE (30) 
k 



6 Alternately, this can be used to fix the Immirzi parameter 7. One can show that 7 = ~~7^- 
Comparing this to the isolated horizons work of Ref. [E3|, if the length scale "lp" in Ref. is 
equal to y/hG, then the result of Eq. (^) agrees with that calculation. 
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is satisfied with the areal radius definition R := |y To leading order in (N), 
the entropy 5, energy (E) and area (A) scale with (N). Therefore we have from 
Eq. @ that (N) < (TV) 3 / 2 , which is satisfied by all (N); this is not a tight bound! 
The result indicates that the ensembles we have considered do not contain black 
holes, which are conjectured to saturate this bound. 

Further refinements of this model, including the implementation of admissibility 
conditions due to the topology of the boundary and gauge invariance, may be found 
inRef. §|. 

5. Discussion 

In this work we have shown that a straightforward application of the ideas of 
statistical mechanics to a gravitational model yield results strikingly similar to black 
hole thermodynamics. In particular, using the grand canonical partition function 
we have shown that in the low temperature, macroscopic limit the system obeys the 
Bekenstein-Hawking entropy relation and the Bekenstein entropy bound. A new 
result is that the system condenses in the lowest level state in the low temperature, 
macroscopic limit. 

Unlike the previous work in spin net quantum gravity, this work begins with the 
quantum gravitational statistical mechanics based on the Boltzmann weight e~ /3H . 
The pathbreaking work in Refs. p7|-|3Cfl based the statistical weight on the concept 
of the "area ensemble," in which the Boltzmann weight is based on the extensive 
parameter area and its conjugate surface pressure. Here, we simply explore the 
statistical mechanics of the Lorentzian gravitational model defined in Section || It 
is a model in which the degrees of freedom of quantum geometry are encoded in 
non-interacting particles with spin. We do not specify black hole or isolated horizon 
boundary conditions and yet observe much the same thermodynamics. 

There are several technical remarks which ought to be mentioned. We remind 
the reader that the quasilocal operator used in this paper is a candidate opera- 
tor. In fact, in Ref. [ pl| there are two inequivalent expressions for the boundary 
Hamiltonian corresponding to two separate normalizations. In one, the operator 
is expressed in terms of the area operator (which provides the quantization of the 
\/|<7afc| in the lapse of Eq. (|l|)). This is the operator used in this paper. However, 
the spectrum of the quasilocal energy operator is modified when vertices lie on the 
bounding surface.^ The spectrum for higher valence vertices (type (ii)) is given by 

p am \ \r Pv(Pv + 2) + n v (n v + 2) - z v (z v + 2) 

EdY,(N) | s) = to 2_ r I s). 

«ed£nG [%Pv(pv + 2) + 2n v (n v + 2) - z v (z v + 2)] 4 

This could be incorporated into the present context by adding another chemical 
potential. Simple statistical arguments suggest that in a fluctuating surface the 
dominant contribution would come from the bivalent intersections used here. 

The second operator in Ref. J3l|] uses a distinct quantization of the classical 
Hamiltonian. The essential difference is that the determinant of the metric is 

7 There is an important caveat to this modified spectrum. In Ref. j^J the quasilocal energy 
operates on states based on networks in the volume and the boundary. After a careful study of 
the quantum boundary conditions and Hilbert space, the Hamiltonian may operate only on states 
based on graphs in the volume. The area operator in Ref. ]jo| has this property. If this is case, 
then the work done in this paper is far more general. It would be valid for all vertices, regardless 
of valence. 
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incorporated using the volume operator. Both operators are well-defined but on 
different spaces. This operator requires that the intersections between the spin 
network state and the boundary are all valence one (type (i)), meaning that graphs 
on which states are based cannot have edges tangent to the surface 9E. This 
spectrum of this operator is 

E d vW) I Si) = mJ2^ nv{n : + 1) I Si), (31) 

v 

where A„ is the eigenvalue of the volume 0, . On account of the state restriction, 
the numerator always has the simple form of the area operator; the operator is 

~ 2 

proportional to A v /V v . The eigenvalues X v can be computed using recoupling 
theory p5| . However, since one can always find a surface tangent to a given edge 
of an embedded graph, the requirement of "no tangent edges" seems too strict to 
apply to general boundaries. 

Throughout the statistical mechanic computation we set the lapse to 1 on the 
bounding surface <9E. Since the lapse is fixed and constant by the boundary con- 
ditions, our choice is one solution. This choice, however, reflects the reparameter- 
ization freedom of time and the selection of which temperature we choose to use 
to investigate the model. Thus, the present statistical mechanics is for an observer 
at the spatial boundary with unit lapse. It would be interesting to investigate a 
model with another boundary and a parameter which measured its size so that one 
could compare solutions as in Ref. Q. 

We close with three comments on the wider implications of this work. In the 
low temperature, macroscopic limit the boundary particles condense into the lowest 
level. Each spin network edge contributes roughly a Planck area (V3£ 2 ) and has 
a degeneracy of 2. In this limit, the model looks like a collection of Planck area 
plaquettes which can each be in one of two states, just like a bit. This is precisely 
Wheeler's "It from Bit" picture of black hole entropy motivated by information 
theory ||. 

In this same limit, variations in fundamental geometry arise from fluctuations 
in the number of spin-1/2 edges. This effectively reproduces the same parameters 
as the Bekenstein-Mukhanov model |l3j - the surface area is a multiple of a funda- 
mental area. This has far-reaching consequences for the spectrum emitted from the 
boundary. Such a configuration would emit a discrete line spectrum p4[ . Thus, if 
edge creation and anhilation is statistically favored, as our analysis suggests, then 
the dominant contribution to radiation would be the Bekenstein-Mukhanov line 
spectrum. 

Finally, it is surprising that such a simple system of independent "bits of geom- 
etry" has the essential elements to reproduce the thermodynamics of the gravita- 
tional field. Many other models and theories have also derived the same thermo- 
dynamics. In string theory, a 1-dimensional gas of D-branes gives the entropy for 
extremal black holes. In Ref. |50| a U(V) Chern-Simons theory gives similar results 
for spherically symmetric black holes. Considering the huge differences in under- 
lying theories, it seems unavoidable that there is a strong element of universality 
of these results. It appears easy to derive the thermodynamics of the gravitational 
field from vastly different theories of fundamental geometry. In this sense it seems 
that matching thermodynamic results is not a stringent test of quantum gravity. 
With an array of theories giving the "correct" thermodynamic limit, we may have to 
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look to new principles to select the correct quantum theory which underlies general 
relativity. 
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